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Abstract 

Magnetic dissipation through decay of Alfven waves in the magnetar magnetosphere 
is discussed. Transport of magnetic fields in the star leads to dissipation of the 
magnetic energy through either direct internal heating or transferring of the energy 
in waves that decay in the magnetar magnetosphere. In the latter case, the Alfven 
waves are excited by crust dislocations or elastic waves underneath the star's surface. 
It is suggested that these Alfven waves can decay into ion sound waves which can be 
effectively damped leading to strong plasma heating. Hot plasmas expand producing 
transient X-rays. 

Key words: pulsars: general: X-rays - magnetars. 



1 Introduction 



Soft gamma-ray repeaters (SGR) and anomalous X-ray pulsars (AXP) are 
believed to be magnetars, wliicli are extremely strongly magnetized, slowly 
rotating neutron stars with a rotation period P = 6 — 8 s and a magnetic 
field well exceedi ng the quantum electrodyna mic (QED) critical field, B ^ 
^ 4.4 X 10^ T flThompson fc Duncan) . flM . Both SGR and AXP emit per- 
sistent thermal X-rays modulated by magnetar's rotation, with a luminosity 
much higher than their spin-down power, suggesting th at the emission is pow- 
ered by the magnetic energy not the rotation energy (IThompson &: Duncan! . 
19951 ). SGRs emit b ursts of soft gamm a-rays (or X-rays) as well as persis- 



tent, pulsed X-rays (IWoods et all . l2006l ). Although the prominent feature of 



AXPs is their persistent X-ra ys, at l east one of them has been detected to 
have SGR-like burst emission (IKaspil . l2004j ) . It is generally thought that radi- 
ation processes in magnetars are predominantly driven by magnetic field decay 
in the neutron star's interior or crustal region, either through direct internal 
heating as the result of a dissipation process or thr ough damping of waves that 
are excited by the release of the magnetic stress (ICumming. Arras. Zweibell . 
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2004J ). This latter mechanism is considered here since no efficient internal dis- 
sipation has been identified so far. An enormous magnetic stress built up in 
the crust as the result of the magnetic transport processes in the star causes 
crustal dislocations or disturbances (e.g. shear waves) that excite fast-mode 
and Alfven waves transferring a significant fraction of the magnetic energy to 
the magnetar magnetosphe re (IBlaes. Blandford. Goldreich. &: Madaul . Il989l : 
Thompson &: Duncan! . Il995l ). Since low-frequency fast-mode waves or Alfven 
waves generated in a superstrongly magnetized plasma are subject to negli- 
gible damping, they mus t unde rgo a cascade through wave-wave interactions. 
Heyl fc HernquistI (119991 . l2005l ) discussed that fast-mode waves may steepen 
and breakdown leading to dissipation of wave energy. We consider here that 
shear waves propagate in a horizontal direction in the crust region and that 
mode coupling produces predominantly Alfven waves. 



Cascades of Alfven waves through three-wave interactions from smaller to 
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This idea was ex- 



tended to the case of a superstrong magnetic field and it was suggested as 
an import ant channel for dissipation of Alfven waves in t he magnetar magne- 
tosph ere (IBlaes. Blandford. Goldreich. fc Madaul . Il989l : iThomspon fc Blaed . 
19981 ). Despite the extensive discussions in both the contexts of the ISM 
and magnetar, the specific dissipation mechanism is not well understood. 
Three-wave interactions in the incompressible MHD approximation lead to 
anisotropic turbulence with the cascade in parallel (to the magnetic field) 



Bhattacharjee &: Ng| . l200ll : lGaltier et al 



wave i iumber suppressed (|Shebalin et all . 119831 : 
2OO2I ; iLithwick fc Goldreichl . |2003| ) , while wave damping through kinetic pro- 
cesses is possible only when a cascade to large parallel wavenumbers occurs. 
Although three-wave interactions involving the fast mode allow cascades in 
both parallel and perpendicular wavenumbers, t hese processes are important 



only below a critical frequency (jChandranl . l2005l ) 



In this paper, decay of Alfven waves through three-wave interactions into 
damping ion sound waves is considered as a possible dissipation mechanism 

in the magnetar magnetosphere. Decay of an Alfven wave into the ion sound 

mode and a backward propagating Alfven wave was suggested by lGaleev fc Oraevskii 
( 1l963l ) and discussed widely in terms of parametric or modulatio r ial instability 
with the applicatio n to solar fiares or the solar wind (jGoldsteinl . Il978l : iDerbyl . 



1978; Miller 



I991I ). In general, when density fiuctuations are allowed, a pump 
wave (not necessary in a natural mode) decays into two side-band plus a den- 
sity modulation. In the low Vs/va -C 1 limit, where Vg is the ion sound speed 
and va is the Alfven speed, one side-band is in the Alfven mode and the density 
fiuctuations are in the ion sound mode. The other side-band is not in a natural 
mode and thus, it is not important. This process reduces to three- wave decay 
( Galeev fc Oraevskii . [l963.) . with the growth rate for the instability (for the 
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side band) being the most effective in the low Vs/va ^ 1 hmit. However, since 
one has a typical Vs/va ~ 0.1 in solar ffares and ~ 0.5 in the solar wind, its 
applicability is limited. Here it is suggested that such process may be efficient 
in magnetars where we have Vs/c ~ 10~^. Since the presence of a propagating 
ion sound wave requires nonequibrilium between electron temperature and ion 
temperature, a condition which is generally n ot satisf i ed for magnetars, the 



sound mode can be treated as a quasimode (IMelrosd . Il986l ) and an Alfven 
wave can decay into an Alfven wave propagating in the opposite direction 
and an ion sound wave that is instantaneously damped. This process leads to 
conversion of the energy in the incoming wave to thermal kinetic energy in the 
plasma. It is proposed here that plasma heating through such process leads to 
X-ray emission. 

In Sec. 2, theory for generation of Alfven waves is outlined and decay of Alfven 
waves through three-wave interactions involving ion sound waves is discussed 
in Sec. 3. A model for thermal X-ray emission due to dissipation of these waves 
in magnetars is discussed in Sec. 4. 



2 Excitation of Alfven waves 



Magnetic transport in the neutron star can lead to magnetic stress being built 
up in the crust which causes fractures or dislocations in the crust. Such events 
generate disturbances propagating, for ex ample as shear waves that can couple 



to th e Alfven wave in the magnetosphere (IBlaes. Blandford. Goldreich. fc Madad . 



19891 ). Although shear waves can couple to both the Alfven mode and fast 
mode, here we discuss only the former where the transmitted energy is as- 
sumed to be mainly in the Alfven mode. One may treat the star's surface 
as a sharp boundary, a reasonable approximation because the scale height of 
the star's atmosphere is about a few cm, much less than the relevant wave- 
length. The typical frequency of shear waves in the crust can be estimated 
from oj ~ Vsh/ld = 10^(K/i/10^ ms~^)(/d/10 m)~^ s~^, where Id is the typical 
length scale of the displacement of crust movement, Vgh = {f^/ pY^"^ is the shear 
velocity, /i ~ 10^^ J m~^ is the shear modulus, and p ~ 10^^ kg m~^ is the crust 
density. 

When the surface is treated as a sharp boundary, Snell's law implies that 
kshsinash = fcsina, where ksh and agh ~ 7r/2 are wave number and incident 
angle (relative to the surface norm) of the shear wave, and k and a are the 
corresponding wave number and incident angle of the transmitted Alfven wave. 
Since the dispersion of the Alfven mode can be written as a; = |fc|||fo, where 
vo = va/ (1 + Va/ c^y^"^ , the propagation angle can be determined from cos 6* = 
k\\/k. In the magnetar's strong magnetic field, one has the limit va c and 
this implies that Vq ~ c. Using u ~ \k\\\c ~ kghVsh, where the phase speed 
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^ Crust 
ksh 

Fig. 1. Coupling of shear waves to the Alfven mode. The surface is treated as a 
sharp boundary so that the Snell's law applies: kgh smash ~ fcsina. Note that 
near the surface the propagation of the Alfven mode is mostly perpendicular, with 




of the shear wave is assumed to be ~ Vgh, one obtains cos 6* ~ K/i/c Thus, 
the transmitted Alfven wave propagates nearly perpendicular to the magnetic 
field at ^ 7r/2 — Vsh/c (as shown in figure [I]). The mode coupling (to Alfven 
waves) is efficient only in the surface area where the magnetic field lines have 
an angle ~ Vsh/c with respect to the surface normal direction. 



A relevant issue here is how the Alfven waves dissipate in the magnetosphere. 
The Landau damping rate is proportional to u^u/Qi)'^, which is negligible 
for the frequency u being much lower than the ion cyclotron frequency 
In the absence of direct damping processes, the Alfven waves must undergo 
cascades through three-wave interactions. A widely discussed process is the 
three-wave interaction involving three Alfven waves (denoted by A, A', and 
A'): A —>■ A' + A" or A + A' ^ A" . Although a cascade in fcy is possible in 
principle in the presence of a natural broadening in the three-wave resonance, 

ly in fc I leading to anisotropic Alfven i c tur- 
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of Alfven and fast-mode (F) waves through A ^ A' + F. Compared to pure 
Alfvenic turbulence, three- wave processes involving the fast mode are more ef- 
ficient for frequenc y below a critical frequency Ur (corresp onding to a critical 
wavelength Vo/wc) (jChandranl . l2005l ; ILuo fc Melrosd . l2006l ). 
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Fig. 2. Kinematic conditions oi A ^ A' + S. Since va S> Vs one has /cy — fcy and 
A;|[ « 2fc||. 

3 Decay into ion sound waves 

Alfven waves can decay into ion sound waves through three wave interactions 
A ^ A' + S. The ion sound mode, like the normal sound wave, is longitudinal 
and is heavily damped unless the ion temperature is much higher than the 
electron temperature. The process is efficient in the low (3s = Vs/vq ^ 1 limit, 
an approximation relevant for the magnetar. Thus, this process provides an 
effective way to damp the Alfven wave. 

3.1 Kinematic conditions 

We consider the weak turbulence approximation in which the wave energy 
and momentum are conserved in each interaction. The energy and momentum 
conservation here, also referred to as the three-wave resonance condition, is 
described hj u = u' + uo" and k = k' + k", where {u, k), {u', k') and (cj", k") 
represent the frequencies and wave vectors of the incoming (A), scattered 
[A') waves and the ion sound {S) wave, respectively. The condition on the 
wave vectors is shown in figure [2J The dispersion of the ion sound wave is 
u" = \k'^\vs with Vs = {KBTjmpY/^ ^ 3x 10^ (Te/lO^ K)i/2 s'S where 
is the electron temperature and we assume ions are protons. In the following 
we assume the incoming Alfven wave propagates forward with respect to the 
magnetic field, fcy > 0. The only permitted processes are those with fcy < 0, 
corresponding to an Alfven wave decays into a sound wave and an oppositely 
directed Alfven wave. From the resonance condition one has 2fc|| = {l+j3s)k" ^ 
k" for Ps ~ Vs/c ^ 1. The decay rate can be determined by the growth rate of 
the backward propagating wave generated from the process. In the following 
we derive the growth rate in the broadband approximation where the relevant 
waves have a broad range of frequencies and wave vectors. 
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3.2 Kinetic equations 



In the broadband approximation, a collection of waves in a particular mode 
can be described by the wave occupation number. Let NA^k), A^^/(fc') and 
Ns{k") be the occupation numbers of the incoming Alfven wave, scattered 
Alfven wave and sound wave. One can write down three kinetic equations 



that determine the rate of change of these three occupation numbers ( iMelrose 



19861 ). To determine the efficiency of the decay we consider the evolution of the 
scattered Alfven wave. Since ion acoustic waves are heavily damped, one can 
neglect A^^. For /?s -C 1, the kinetic equation for the scattered wave {A') due 
to the processes A A' + S and A' ^ A + S can be written approximately 

as 



d\nNA'{k') f dk 



dt 



UJ 



X 



_d_ 

dkw 



UAA's{k,k\k")NA{k) 



where fcy ~ —k'^\, uaa's = f>'Rs\^i^'j^"*iC(iji\'^ / i^f,^^'^") is the probability for the 
transition A ^ A' + S, aiji is the quadratic response tensor, Rs is the ratio of 
the electric energy to the total energy in the sound wave, and the correspond- 
ing ratio for the incoming and scattered waves is taken to be Ra = Ra' = 1/2- 
The three vectors e, e' and e" are the polarizations of the three waves con- 
cerned. The evolution of the scattered wave depends only the occupation num- 
ber of the incoming wave. 



3.3 Growth rate 



The growth of the backward propagating Alfven waves can be calculated from 
([1]). The three- wave probability can be derived assuming that electron mo- 
tion in two high frequency waves (cu, u') is treated in the cold plasma ap- 
proximation. Specifically, we assume the following conditions are satisfied: 
ui" /k'^^ = -C f e -C uj/k\\ = Vq, where Ve is the electron th ermal speed. T he 
quadratic response tensor has the following approximation (IMelrosd . 119861 ) 



a,,s^^-^r,A^)k':x'^^\uj",k"), (2) 
Tij^bibj + i—eijsbs, (3) 
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where x^^^'' ~ ^pe/^"'^'^e electronic contribution to the longitudinal 

component of the dielectric tensor, b = B/B is the magnetic field direction, 
ujpe is the electron plasma frequency, Sijg is the permutation symbol (1 for 
even permutation, —1 for odd permutation and otherwise), and Qe is the 
electron cyclotron frequency. The strong magnetic field approximation implies 
that only terms 0123 = —0^213 are important. Since the component ([2]) corre- 
sponds to the order l/fig terms in the quadratic response tensor, electrons 
and positrons do not cancel in three-wave interactions. In deriving the three- 
wave probability the Alfven waves are assumed to be linearly polarized with 
e = (cos 0, sin0, 0) and e' = (cos 0', sin 0', 0). The wave vectors are written 
as A; = (A;_|_ cos 0, A;_L sin 0, fc||) and k' = (fc^ cos0', fc^ sin0', In the strong 
magnetic field, one has Rs ^ A;"\2/(2cj2.) with LUpi the ion plasma frequency. 
Using ([2]), one obtains 



xsin2(0'-0), (4) 

where Ve ~ 2.8 x 10^^^ m is the classical electron radius. Here we assume that 
the relevant waves are linearly polarized. For the three-wave probability to be 
nonzero, the two linear polarizations must be noncoplanar (0' 7^ 0). Denoting 
the left-hand side of ([1]) by F^/ one obtains 



d 



X / dkt k\ 



dh 



(5) 



where N^ik) = NA{k^,k\\) is assumed to be axially symmetric. The growth 
rate is sensitive to the spectrum of the incoming Alfven waves, a result simi- 
lar to that derived by I Cohen &: Dewarl (119741 ) in the MHD approximation in 
the special case where the Alfven waves are circularly polarized propagating 
parallel or antiparallel to the magnetic field. 

The propagation of Alfven waves emitted from the surface due to coupling of 
the shear waves is oblique k± ^ fcy (Sec. 2), implying Na is anisotropic in the 
k space. As an example, we consider 



NAik±,k\\) oc ku-^S{k±-rikii), 



(6) 



where the parameter t] characterizes the anisotropy, which can be taken as ~ 
c/Vsh ^ 1 for the emission mechanism discussed in Sec. 2. This corresponds to 



7 




Fig. 3. Plasma heating due to decay of Alfven waves near the surface. The sohd 
arrows indicate the group velocity of the relevant waves, which are directed either 
parallel or antiparallel to the magnetic field direction. The backward propagating 
wave may be reflected from the surface. So, the waves may be trapped near the star's 
surface. The ion sound wave is assumed to be absorbed instantaneously leading 
to plasma heating. The shade represents a hot plasma (proton+electron) cloud 
emitting X-rays. 

an energy spectrum Ea{uj) ~ ci;" of the incoming wave. It can be shown that 
the positive growth (corresponding to decay of the incoming waves) requires a 
spectrum with a positive slope. This condition can be understood as follows. 
For a transition A A' + S, there is an inverse process A' ^ A + S that 
removes the A' waves. For A ^ A' + S to he dominant over A' ^ A + S, 
an inversion (i.e. a positive slope) in the occupation number Na is required. 
Assuming that the spectrum is peaked at a; ~ cuq with the positive slope 
characterized by a power index a > 0, the growth rate is derived as 

_ a{a + 1) fujpi^Y ( ^Ba V' , . 

uj' ATi^ \uoJ \ B ) ' ^ ' 



Although the parameter (3s is not present explicitly in the growth rate, the 
condition (3,^ is needed so th at the three-wave resonance condition can be 
satisfied (ICohen &: Dewarl . Il974l ) . 



4 X-ray emission from plasma heating 



In the following we specifically consider SGRs, which are characterized by 
sporadic outbursts as well as occasional giant bursts. The burst emission from 
SGRs is generally considered to be magnetospheric origin, while the persis- 
tent, pulsed emission may be attributed to surface emission. Normal outbursts 
typically last weeks or months and consist of many ordinary, repetitive short 
bursts that last about 0.1 s, which sets a lower limit to the damping time for 
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the incoming Alfven wave. This damping time can be determined from the 
growth time of the backward propagating wave (cf. Eq. [7]). 



To estimate the growth rate one needs to assume a specific form of the in- 
coming wave spectrum which should be determined by the shear waves in the 
crust. So far, no detailed model is available for the calculation of the latter. 
However, one may make an analogy with the energy spectra of shear waves 
in earthquake s by assuming a spectrum peaked at a characteristic frequency 



given by uq ( iBlaes. Blandford. Goldreich. fc Madaul . Il989l ). The decay pro- 



cess discussed in Sec 3 occurs for uj' < uq. We assume that emission of Alfven 
waves from the surface is accompanied by injection of plasmas of protons 
plus electrons (or/and positrons) into the magnetosphere. The plasma den- 
sity is not well constrained but its minimum density should be about the GJ 
density ugj = ^nsoB/eP ^ (1.4 x 10^^ m-3)(P/5 s)-i(5/10^^ T), where we 
ignore the cosine angle between the magnetic field and the rotation axis. As- 
suming Up ^ ncj ~ 10^^ m~^, k'^ = 300fc||. about the GJ density, a = 2, 
and SBa/B = 5 x 10~^, the growth time is estimated to be 1/Ta' ~ lO^'^s 
for o^o ~ 10^ s~ ^. This time is much shorter than the cascade time through 



A ^ A' + A" (ILuo fc Melrosd . 120061 ) . implying that the decay occurs at a 



much faster rate than pure Alfvenic cascades. The backward propagating wave 
generated from the process may be reflected from the surface and propagates 
forward undergoing further decay. As a result, the waves are trapped close to 
the surface (at a distance As = c/Ta' ~ 2 x 10^ m from the surface) without 
propagating to the opposite foot points (as shown in figure [3]). The trapped 
waves continue to decay into ion sound waves which are absorbed by ther- 
mal electrons. As a result, the plasma can be rapidly heated leading to X-ray 
emission. Since the fraction of the wave energy dissipated in each interaction 
is u!"/u! ~ Vs/c ~ 10~^, the time for the trapped waves to be absorbed is 
{c/vs)/Ta' ~ 0.1 s, which is comparable to the typical time scale of individual 
bursts. 

An attractive aspect of the mechanism discussed here is that the X-ray emis- 
sion from the decay of Alfven waves is intrinsically variable on the time 
1/Ta'- Such time scale may b e relevant for variabilities observed in giant flares 



( IWatts fc Strohmayen . l2006l ). However, to assess its relevance a quantitative 
model for the emission region is needed, which is beyond the scope of this 
paper and will be discussed elsewhere. 



5 Conclusions and discussion 



Transport of magnetic fields in the neutron star causes magnetic dissipation 
both in the core and outside the star. The former leads to heating which may 
be radiated through neutrino emission or partially conducts to the surface 
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producing thermal X-ray emission. In the latter case, the magnetic energy 
is efficiently transported to the magnetar's magnetosphere in Alfven waves. 
These waves decay through three-wave interactions. It is suggested here that 
because of the low /^^ -C 1, the process involving a Alfven wave decaying into 
ion sound waves that are then damped is the most efficient way to dissipate 
the magnetic energy. In each interaction, the Alfven wave cascades toward 
the small wavenumber regime. Since the process is intrinsically dissipative 
converting the energy in the wave to kinetic energy in the plasma, such 'in- 
verse cascade' can be understood in analogy with the Langmuir condensation 
due to nonlinear Landau damping which tends to pump the Lanmguir wave 
from larger to smaller wavenumber (IMelrosd . Il986l ) . The growth rate for the 
scattered Alfven wave, which quantifies the efficiency of the process, is calcu- 
lated in the random phase approximation that is applicable provided that the 
waves generated have a broad range of frequencies and wave numbers. This 
approximation is relevant since the crust dislocations or fractures may produce 
many shear waves that are not related in phase (i.e. incoherent waves). It is 
shown here that the absorption time can be comparable to or shorter than the 
typical duration of individual outbursts, suggesting that such process can be 
important for transient X-ray emission. 



Note that the decay process considered here was widely discussed in the MHD 
approximation in the context of modulational instability or parametric de- 
cay with the ap p licati o n to the sola. r wind or solar flares (iLashmore-Davied . 



19761 : iGoldsteinl . Il978l : iDerbyl . Il978l ). There were also recent discuss i ons o f 



this process using numerical simulation (jPel Zanna. Velli. Sz Londrilld . l200ll ). 
These discussions mostly concentrated on circularly polarized Alfven waves 
and assumed that the incoming wave is monochromatic with the growth 
rate estimated to be Ta'/coo ~ [5 B a / B) / jSlf " ^ in the limit (3s = Vs/va <^ 1 



( iGaleev fc Oraevskiil . 1 19631 : IGoldsteinl . 1 19781 : iDerbyl . Il978l ). It is debatable 



whether such result is relevant for the solar wind and solar flares because 
in both cases, the incoming waves should be broad band and it may not be 
appropriate to treat the incoming waves as a single coherent wave. Further- 
more, in both cases the parameter (3s = Vs/va is moderate. At a relatively 
large (is the waves concerned (the relevant side-band waves) are away from 
the corresponding natural modes and as a result, the growth rate is consider- 
ably reduced. It is shown here that the decay process is favored in magnetars 
and can be an effective mechanism for magnetic dissipation required for the 
interpretation for the observed transient X-ray emission. 
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